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Narrative 


Philip  Holmes  is  continuing  to  study  the  mathematical  foundations  of  our  general 
approach  to  the  turbulent  boundary  layer:  projection  of  the  Navier-Su^es  equations  onto 
(relatively)  low  dimensional  subspaces  spanned  by  empirical  eigenfunctions.  In  particular, 
in  a  brief  paper  prepared  for  the  ICM  *90  Proceedings  (Kyoto,  Japan),  Berkooz,  Hcrimes 
and  Lumley  show  that,  on  the  average,  the  coupling  terms  between  resolved  modes 
streamwise  variadon  and  unresolved  nodes  camiot  transfer  energy  from  the  latter  to  the 
former,  even  if  relatively  low  wavenumbers  are  among  the  laner,  neglected  mndes. 
Holmes  also  continued  his  study  of  symmetry-induced  heteroclinic  cycles,  their  relation  to 
the  bursting  process  (and  more  generally  to  intertnittency  in  turbulence  production)  and 
their  behavior  in  the  presence  of  multiple  local  and  global  instabilities  in  phase  space. 

Peter  Swan  is  currently  investigating  (with  P.  Holmes),  the  large-time  behavior  of 
systems  that  can  develop  fine  structure.  In  particular  they  are  looking  at  some  visco-elastic 
models  relevant  to  structural  phase  transitions  -  for  these  investigations  they  have  just  been 
awarded  supercomputer  resources  by  the  NSF  supercomputer  -  allocations  committee. 

Swan  received  an  MSI  Graduate  Fellowship  for  the  academic  year  1990-91. 

During  the  past  nine  months  Mark  Myers  has  primarily  provided  technical  support 
for  the  development  of  the  kaos  dynamical  system  toolkit,  developed  by  Drs. 
Guckenheimer  and  Kim.  During  this  period  several  parts  of  the  program  were  enhanced,  a 
delivery  version  of  the  package  was  generated  and  a  users’  guide  was  prepared.  The 
addition  of  new  code  to  the  program  in  the  fall  was  substantial;  for  example,  a  complete 
conversion  of  the  eigenvalue  and  eigenvector  computation  routines  in  the  EISPACK  library 
from  FORTRAN  into  the  C  language  was  performed  so  that  these  well-established 
algorithms  could  be  included  in  kaos.  Since  the  kaos  package  was  made  publicly 
available  in  February,  a  number  of  error  reports  have  been  received  and  appropriate  action 
taken  by  various  members  of  the  kaos  suppon  group.  It  is  expected  that  continued 
suppon  of  this  kind  for  the  kaos  package  will  be  required  in  the  future. 

Concurrent  with  the  kaos  development  work,  the  Applied  Mathematics  computing 
center  facilities  were  greatly  enhanced  and  suppon  by  the  research  grant  allowed  Myers  to 
participate  in  this  expansion.  The  center  currently  suppons  seven  Sun  workstations 
available  for  use  by  faculty  and  students,  two  of  which  were  purchased  by  our  group. 
These  computers  are  used  frequently  by  members  of  the  turbulence  research  group;  system 


support  and  maintenance  is  crucial  if  we  are  to  insure  that  the  group  continues  to  benefit 
from  these  resources. 

Sidney  Leibovich  was  elected  Chairman  of  the  U.S.  National  Committee  on 
Theoretical  and  Applied  Mechanics  (A  committee  of  the  NRC  and  Nat  Academies).  1990- 
92.  He  was  also  elected  U.S.  Delegate  to  the  lUTAM. 

Sue  Campbell  (in  collaboration  with  Phil  Holmes)  is  just  finishing  up  work  on  the 
interactions  between  structurally  stable  heteroclinic  cycles  and  travelling  waves  in  a  three 
mode  truncation  (modes  1, 2  and  4)  of  the  model  studied  by  Aubry  et  aL;  they  have  a  paper 
ready  for  submission  entided  ''Bifuicarion  firom  symmetric  heteroclinic  cycles  with  three 
interacting  modes". 

The  purpose  of  this  work  was  to  better  understand  the  heteroclinic  behavior  seen  in 
the  many-mode  models  by  Aubry  et  al.,  which  is  thought  to  be  related  to  the  bursting 
phenomenon  observed  in  turbulence.  As  the  heteroclinic  cycles  often  appeared  between 
fixed  points  in  the  2:4  subspace  of  the  larger  systems,  it  was  apparent  that  in  order  to 
understand  the  behavior  in  the  larger  systems  one  must  first  understand  the  dynamics  in 
this  subspace. 

The  1 :2:4  mode  system  was  chosen  to  be  studied  in  the  hope  that  it  would  be  large 
enough  to  retain  much  of  the  behavior  in  the  larger  systems,  but  small  enough  to  be  tackled 
analytically.  The  former  was  certainly  true,  for  much  rich  behavior  reminiscent  of  that  seen 
in  the  higher-mode  models  was  observed  through  numerical  experiments.  The  latter  was 
somewhat  true,  and  with  the  aid  of  MACSYMA  some  results  on  the  existence  of  the 
heteroclinic  cycles  and  their  bifurcation  to  modulated  travelling  waves  were  established. 

Andrew  Poje  is  currently  working  (with  John  Lumley)  on  the  eduction  of  coherent 
structures  in  wall  bounded  flows  by  analytical  methods.  This  work  is  taking  place  m  a  two 
pronged  fashion.  The  problem  as  originally  formulated  by  Lumley  (1969)  has  t«en  solved 
numerically  for  both  channel  and  open  boundary  layer  flow.  An  extensive  data  set  is 
available  for  the  calculation  of  coefficients  in  the  dynamical  systems  derived  by  Aubry  et  aL 
(1987).  A  study  of  the  dynamics  of  systems  formed  using  analytic  eigenfunctions,  and  a 
comparison  with  and  extensions  to  the  work  of  Aubry  et  al.,  is  currently  taking  place. 
Along  these  lines,  the  recent  acquisition  of  a  large  set  of  Proper  Orthogonal  Decomposition 
eigenstructures  from  DNS  data  of  Kim  and  Moin  allows  for  direa  comparison  of  PC®  and 
linear  energy  eduction  techniques.  Preliminary  results  of  the  scaling  behavior  of  the  energy 


content  of  structures  derived  by  both  methods  is  available  in  a  note  being  composed  by 
Berkooz  and  Poje  (see  below). 

Poje  is  also  pursuing  another  aspect  of  linear  energy-method  stability  analysis;  this 
focuses  on  extensions  of  the  technique  to  more  complicated  flows.  A  simple  naodel  for  the 
effects  of  riblets  on  the  growth  and  form  of  eigenstnictures  is  under  investigation.  This 
model  is  a  straight-forward  extension  of  previous  work  to  flows  with  mean  profiles  that 
depend  on  both  wall-normal  and  spanwise  coordinates.  A  more  complicated  extension  ot 
the  eduction  method  to  the  case  of  the  thermal  boundary  layer  is  under  investigation.  This 
flow,  with  applications  to  boundary  layer  control,  is  of  acute  meteorological  interest  U  is 
hoped  that  the  relatively  simple  eduction  method  under  consideration  will  allow  the 
application  of  dynamical  system  techniques  to  models  of  the  atmospheric  boundary  layer. 

Scott  Jones  (in  collaboration  with  Leibovich)  is  currently  investigating  shear 
dispersion  by  chaotic  advection.  He  has  performed  preliminary  calculations  of  the  sdiring 
of  fluid  particles  in  a  flow  composed  of  a  superposition  of  a  second  order  irrotational  wave 
solution  plus  a  linear  shear.  He  has  demonstrated  that  this  flow  can  generate  chaodc 
particle  trajectories.  Starting  with  particles  initially  configured  as  a  circle,  the  flow  given  by 
a  two-mode  approximation  plus  shear  causes  exponential  growth  in  a  measure  of  particle 
separation. 

Jones  and  Leibovich  are  also  applying  the  ideas  from  Langmuir  circulation  to  the 
wall  region  of  a  turbulent  channel  flow.  This  is  done  using  the  Generalized  Lagrangian 
Mean  formulation.  Lumley  and  Leibovich  have  concluded  that  the  wall  region  of  a 
turbulent  channel  flow  or  boundary  layer  is  very  similar  to  the  ocean  surface  mixed  layer, 
but  upside  down.  The  turbulent  fluctuations  in  the  wall  region  generate  a  pseudo¬ 
momentum  (related  distandy  to  the  Stokes  drift)  just  as  the  surface  gravity  wave  motions  do 
in  the  ocean  surface  mixed  layer.  We  expect  the  same  son  of  instability,  by  the  same 
mechanism  (stretching  of  perturbed  cross-stream  vorticity  by  shear  in  the  Stokes  drift), 
leading  to  longitudinal  rolls  in  both  cases.  We  identify  these  rolls  as  the  coherent  structures 
observed  in  the  wall  region,  and  implicated  in  the  bursting  process.  Lumley  has  generated 
approximate  forms  for  the  pseudo-momentum  in  the  wall  region,  using  data  from  the 
measurements  of  Favre  et  al .  These  forms  will  be  used  in  calculations  of  the  genesis  of  die 
rolls. 

Stephen  M.Cox  is  also  working  on  Langmuir  circulations,  but  in  a  stratified  ocean, 
with  Sidney  Leibovich.  They  are  lot^ng  at  the  possibility  of  modeling  the  slow  and  long- 


wavelength  modulations  of  the  amplitudes  of  the  rolls  by  using  coupled  Ginzburg-Landau- 
type  equations.  This  son  of  analysis  can  also  be  applied  to  the  slow  and  long-wavelength 
modulations  of  the  amplitude  of  rolls  in  the  wall  region  of  the  turbulent  boundary  layer. 

Cox  is  also  collaborating  with  John  Guckenheimer  on  the  inieracdon  of 
parametrically  excited  surface  waves  in  a  nearly  square  container  subject  to  vertical 
oscillations.  They  are  trying  to  make  sense  of  apparent  disagreements  between  weakly 
non-linear  theory  and  experimental  evidence.  The  equations  describing  this  system  have 
many  points  of  similarity  with  the  equations  describing  the  wall  region  of  a  turbulent  flow, 
but  are  in  some  respects  more  tractable.  It  is  expected  that  insight  gained  in  one  system  will 
carry  over  to  the  other. 

Gal  Berkooz  has  given  primary  attention  during  this  report  period  to  understanding 
and  explaining  the  behavior  exhibited  by  the  low  dimenskmal  dynamical  systems  «»»«***■*  in 
Aubry  et  al.,  in  particular  in  light  of  Moffat's  observations  (in  Whither  Turbulence?). 
Moffat  pointed  out  that  streamwise  rolls  without  streamwise  variation  decay  according  to 
the  exact  equations,  whereas  in  our  model  system  they  are  able  to  maintain  themselves.  To 
this  end  Berkooz  did  the  following: 

He  showed  that  the  notion  of  coupled  and  uncoupled  systems  of  equations  is  not 
determined  by  the  form  of  the  equations  alone,  but  is  also  dependent  on  the  nmm  defined 
on  the  space  in  which  one  is  working.  This  is  the  "mathematical"  observation  which 
answers  Moffat's  criticism. 

He  collaborated  with  J.L.  Lumley  and  P.Holmes  in  showing  that  the  bursting 
behavior  displayed  by  the  low  dimensional  model  system  is  not  an  artifact  of  the  coupling 
introduced  by  the  modeling,  but  persists  in  an  uncoupled  (decaying)  system  They  also 
showed  that  the  production  of  energy  in  these  systems  is  within  the  experimentally 
observed  values  for  the  modes  used  in  the  truncation. 

Berkooz  also  showed  that  several  assumpdons  made  on  an  intuitive  basis  in  the 
work  of  Aubry  et  al.  may  be  justified  formally.  Namely: 

That  the  Heisenberg  model  used  gives  the  correct  dissipation  to  within  a  factor  of 
order  unity,  as  assumed. 

That  the  Leonard  stresses  may  be  neglected  in  the  case  of  modeling  with  no 
streamwise  variation,  as  assumed. 


That  the  previous  result  holds  for  an  arbitrary  number  of  eigenfunctions  when  no 
streamwise  variations  arc  present 

Berkooz  gave  a  physical  meaning  to  the  notion  of  models  with  no  streamwise 
variation  and  showed  that  in  effect  they  average  the  streamwise  dynamics,  as  conjectured 
by  Holmes. 

Berkooz  also  applied  dynamical  systems  models  to  the  control  of  the  turbulent 
boundary  layer.  As  observed  by  Lumley,  the  dynamical  systems  models  are  the  only 
mechanistic  description  of  the  dynamics  in  the  wall  region  of  the  turbulent  boundary  la^, 
and  as  such  are  the  only  cartdidates  for  establishing  a  smart  control  scheme.  Early  work  of 
Marsden  and  Bloch  showed  that  systems  with  homoclinic  attractors  are  in  principal 
controllable,  with  a  certain  type  of  control  input  Berkooz's  efforts  in  this  direction  were 
concentrated  on  determining  the  feasibility  of  the  control  and  trying  to  understand  the 
possible  gains  in  terms  of  drag  reduction  or  mixing  enhancement  To  this  end  he  did  the 
following: 

He  introduced  the  notion  of  short  term  tracking  time  (denote  T_s).  This  is  a  time 
scale  characteristic  of  the  time  over  which  the  dynamical  system  model  tracks  the  true 
dynamics  accurately.  T_s  is  of  fundamental  importance  in  the  control  application;  it  must 
be  of  the  order  of  the  wall  region  time  scales  to  enable  control. 

Berkooz  showed  that  dynamical  systems  based  on  the  Proper  Orthogonal 
Decomposition  have,  on  the  average,  the  best  T_s  for  a  given  number  of  modes. 

To  determine  the  effect  of  control  actuators  on  the  structure  of  the  turbnlem 
boundary  layer,  Berkooz  suggested  a  numerical  scheme  for  a  direct  pseudo-spectral 
numerical  simulation  .  Work  is  currently  in  progress  in  implementing  this  scheme,  with  D. 
Bond. 


Berkooz  also  submitted  a  proposal  to  the  Center  for  Turbulence  Research  summer 
session  for  determining  T_s  based  on  their  D.N.S.  data.  The  proposal  was  accepted  and 
hopefully  the  work  will  take  place  in  the  summer  of  "90. 

Berkooz  also  generalized  previous  wOTk,  turning  the  dynamical  systems  approach 
to  a  quantitative  predictor,  among  other  pertinent  topics.  To  this  end  he  did: 
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He  generalized  the  work  of  Annbruster.  Guckenheimer  and  Holmes  concerning  the 
structural  stability  of  heteroclinic  connections  in  systems  with  0(2)  symmetry  to  three 
complex  dimensions,  supporting  the  conjecture  that  the  previously  established  results  are 
generalizable  to  higher  dimensions. 

He  recognized  that  to  get  quantitative  predictions  will  require  generating  higher  and 
higher  dimensional  dynamical  systems,  and  that  the  best  way  to  go  about  this  is  to  uiilue 
computer  algebra.  To  this  end  he  initiated  an  interdisciplinary  effort  to  meet  this 
requirement  Participants  include  Professtvs  R.  Zippel  from  Computer  Science,  S.  Rapkin 
of  Computer  Science  and  Y.Y.  Du  of  Applied  Math. 

He  collaborated  with  A.  Poje  in  studying  the  asymptotic  behaviour  of  empirical 
eigenvalues  and  energy  stability  eigenvalues.  This  is  required  if  one  wants  to  assign 
energy  content  to  stability  analysis  eigenfunctions  in  an  internally  consistent  manner.  This, 
in  turn,  is  essential  for  utilizing  stability  analysis  eigenfunctions  instead  of  empirical 
eigenfunctions  in  dynamical  systems  models. 

The  1982  experiment  by  Hopfinger,  Browand  and  Gagne  (HBG;  JFM  125:  505* 
534)  displayed  remarkable  self-organization  of  turbulent  flows  in  a  rotating  tank.  The 
turbulent  flows  form  well-ordered  panems  consisting  of  very  strong  concentrated  vortices 
running  the  length  of  the  tank,  parallel  to  the  angular  velocity  vector.  To  understand  this 
pattern  formation  process,  which  has  implications  for  all  flows  in  which  vOTticity  of  several 
scales  is  present,  Guang  Yang  (in  collaboration  with  S.  Leibovich)  is  planning  a  direct 
numerical  simulation  of  the  flow  between  two  infinite  parallel  planes  with  boundary  foRing 
at  one  plane.  The  Navier-Stokes  equations  are  solved  using  a  pseudo- spectral  method. 

The  code  for  the  numerical  simulation  has  been  written.  Its  accuracy  has  been 
checked  by  computing  the  evolution  of  small  and  finite-amplitude  two-dimensional 
disturbances  of  plane  Poiseuille  flow  and  comparing  the  results  with  their  predicted 
behavior  according  to  the  Orr-Sommerfeld  equation  and  the  numerical  simulation  results  of 
Orszag  and  Kell  {JFM  96:  159-205).  Now  Yang  and  Leibovich  will  simulate  the  HBG 
experiment,  but  at  a  Reynolds  number  lower  than  that  in  the  experiment. 

In  the  course  of  the  1989-1990  academic  year,  Jerrold  Marsden  decided  not  to 
accept  a  permanent  position  at  Cornell,  but  to  return  to  UC  Berkeley,  from  which  he  was 
on  Leave  of  Absence.  We  have  not  yet  decided  how  to  restructure  the  research  program. 
Marsden  has  already  made  a  definitive  contribution,  and  we  would  probably  be  able  to 
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proceed  without  further  input  in  his  area.  Wc  couid  continue  to  work  with  him  at  Berkeley, 
or  we  could  take  on  someone  else  in  the  control  area,  cither  at  Cornell,  or  elsewhere.  We 
are  seriously  considering  taking  on  Anthony  Bloch  (who  did  the  earlier  control  work  under 
Marsden's  direction)  as  a  consultant.  Bloch  is  now  at  Ohio  state.  Electronic 
communication  is  now  such  that  this  would  not  present  a  problem. 
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Research  synopsis 


Current  research  effons  of  Pojc  and  Lumicy  arc  focused  on  accurately  (adequately) 
modelling  the  effects  of  modulated  Reynolds  stresses  on  the  form  of  the  coherent  stnictures 
predicted  by  an  energy  maximization  principle.  We  have  found  that  a  simple  eddy- 
viscosity  relation  does  not  adequately  model  the  stresses  which  may  act  as  productiop  term 
for  the  coherent  structures.  An  anisotropic  eddy-viscosi^  u.nsor  has  bwn  derived  from 
second  order  equations  by  assuming  the  coherent  velocity  field  to  be  an  order  epsilon 
perturbation  to  the  mean  field.  Initial  results  using  this  model  show  coherent  structures 
which  are  quite  similar  to  those  predicted  by  the  Proper  Orthogonal  Decomposition. 

J.  Elezgaray  gave  a  mini-course  on  wavelets  during  the  fall  semester.  The  course 
was  intended  for  people  interested  in  applications  such  as  the  turbulent  boundary  layer. 
The  course  consisted  of  five  lectures,  for  which  lecture  notes  are  available. 

So  far  our  decompositions  have  used  Fourier  modes  in  the  homogeneous  spanwise 
and  streamwise  directions.  These  distributed  functions,  while  energy-optimal  on  average, 
are  not  well  adapted  to  reproduce  spatially  localized  events  such  as  bursts.  In  addition,  one 
of  the  points  on  which  the  work  of  Aubry  et  al  was  criticised,  was  that  the  stnicturea 
produced  were  not  localized  in  space  like  the  structures  observed  in  the  turbulent  bomdrey 
layer.  Wavelet  bases,  having  compact  support,  seem  much  more  appropriate.  J.  Ele^any 
is  currently  studying,  jointly  with  G.  Berkooz  and  P.  Holmes,  applications  of  wavelets  to 
the  study  of  coherent  structures  arising  from  the  Kuramoto-Sivashinsky  equation.  This 
study  is  preliminary  to  a  more  involved  study  of  the  more  complex  rurbulent  boundary 
layer,  permitting  us  to  gain  experience  with  the  applicadon  of  wavelets  in  chaodc  systems. 

Berkooz,  in  collaboration  with  Holmes  and  Lumley,  gave  substantial  attention  to 
understanding  and  explaining  the  behavior  exhibited  by  the  low  dimensional  dynamical 
systems  studied  in  Aubry  et  al.  In  particular,  they  showed  that  several  assurapdons  made 
on  an  intuitive  basis  in  the  work  of  Aubry  etal  may  be  justified  formally,  namely;  that  the 
Heisenberg  model  used  gives  the  correct  dissipadon  within  a  constant  of  order  unity,  as 
assumed,  that  the  Leonard  stresses  may  be  neglected  in  the  case  of  modeling  with  no 
streamwise  variadon,  as  assumed;  that  the  previous  result  holds  for  an  arbitrary  number  of 
eigenfunctions  when  no  streamwise  variadon  is  present;  that  models  with  no  streamwise 
variadon  in  effect  average  the  streamwise  dynamics,  as  conjectured  by  Holmes. 

Berkooz,  Lumley  and  Holmes  submitted  a  proposal  to  the  Center  For  Turbulence 
Research  summer  session  for  detennining  the  short  term  tracking  dme  (T_s)  from  the 
CTR's  data  from  direct  numerical  simuladon  for  channel  flow.  The  proposal  was  accepted, 
and  the  following  work  took  place  during  the  1990  summer  session;  the  DNS  data  sets  for 
the  "minimal  flow  unit  of  wall  turbulence",  studied  by  Jiminez  and  Moin,  were  examined  in 
detail.  The  minimal  flow  unit  is  the  smallest  computational  domain  found  to  sustain  a 
velocity  field  similar  (in  some  statistical  charaaerisdes)  to  a  turbulent  field.  It  was  dweiried 
to  use  this  data  set  since  it  shares  some  key  features  with  the  dynamical  systems  models;  it 
was  also  considerably  simpler  to  access  than  the  full  channel  simulations.  The  time  traces 
were  produced,  for  what  would  correspond  to  the  eigenfunction  coefficients.  These  were 
produced  with  the  idea  of  measuring  T_s,  but  this  proved  to  be  too  involved  for  the  time 
available.  .A  qualitative  numerical  evaluation  of  the  predictions  based  on  the  dynamical 
systems  approach  was  carried  out.  The  main  contribution  of  this  work  was  to  establish  a 
methodology  for  the  evaluation  of  dynamical  systems,  and  to  introduce  some  new 
concepts.  The  amount  of  data  was  insufficient  to  permit  decisive  conclusions. 

Lumley  has  pointed  out  that  the  dynamical  systems  models  are  the  only  mechanistic 
description  of  the  dynamics  of  the  wall  region  of  die  turbulent  boundary  layer,  and  thus  are 
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the  only  candidates  for  establishing  a  smart  control  scheme.  Early  work  of  Marsden  and 
Bloch  showed  that  systems  with  homociinic  attractors  arc  in  principle  controllable,  with  a 
certain  type  of  control  input.  Our  efforts  in  this  direction  are  concentrated  on  determining 
the  feasability  of  control,  and  trying  to  understand  the  possible  gains  in  terms  of  drag 
reduction  and  mixing  enhancement  To  this  end  Berkooz,  Holmes  and  Lumley  introduced 
the  notion  of  short  term  tracking  time  T_s.  This  is  a  measure  of  the  time  over  which  a 
dynamical  systems  model  tracks  the  true  dynamics  accurately.  T_s  is  of  fundamental 
imponance  in  the  control  application,  and  it  must  be  of  the  order  of  the  wall-region  rinM> 
scales  to  make  control  possible.  They  then  showed  that  dynamical  systems  bas^  on  the 
Proper  Onhogonal  Decomposition  have,  on  the  average,  the  best  T_s  for  a  given  number  of 
modes. 


To  determine  the  effect  of  control  actuators  on  the  structure  of  the  turbulent 
boundary  layer,  Berkooz,  Holmes  and  Lumley  suggested  a  numerical  scheme  for  a  direct 
(pseudospectral)  numerical  simulation  scheme.  This  year  David  Bond  and  Berkooz 
the  DNS  channel  code  of  Moin,  Kim  &  Moser  from  CTR,  and  translated  it  into  FORTRAN 
so  that  It  could  be  run  locally.  An  algorithm  was  formulated  for  implementation  of  a 
linearized  boundary  condition  simulating  the  flow  over  a  control  actuator.  The  group 
applied  for,  and  received,  200m  Cray  hours  on  the  NAS  system  (at  NASA/Ames)  for  the 
implementation  of  the  numerical  schone. 

Berkooz  generalized  the  work  of  Armbruster,  Guckenheimer  and  Holmes  regarding 
the  structural  stability  of  heteroclinic  connections  in  systems  with  0(2)  symmetry  to  three 
complex  dimensions,  supporting  the  conjecture  that  the  previously  established  results  are 
generalizable  to  higher  dimensions.  This  is  related  to  the  work  of  (Campbell  and  Holmes 
(see  below). 

Berkooz  recognized  that  to  get  quantitative  predictions  will  require  generating 
higher  and  higher  dimensional  dynamical  systems,  and  that  the  best  way  to  go  about  this  is 
to  utilize  computer  algebra.  To  this  end,  Berkooz  initiated  an  inte^isciplinary  effort 
involving  R.  Zippel  and  S.  Rapkin  from  Computer  Sciences. 

Berkooz  collaborated  with  A.  Poje  in  studying  the  asymptotic  behavior  of  empirical 
eigenvalues  and  energy  stability  eigenvalues.  This  is  required  if  one  wants  to  assign 
energy  content  to  stability  analysis  eigenfunctions  in  an  internally  consistent  manner.  Thii, 
in  turn,  is  essential  for  utilizing  stability  analysis  eigenfunctions  instead  of  empirical 
eigenfunctions  in  dynamical  systems  models. 

Berkooz  studied  the  mathematical  properties  of  the  proper  orthogonal 
decomposition  and  the  phenomena  of  localiz^  coherent  structures  in  homogeneous 
(statistically  translation  invariant)  systems.  He  established  some  results  of  particular 
interest  in  our  context,  regarding  mathematical  properties  of  the  POD. 

Guckenheimer's  work  is  continuing  on  computational  tools  fw  expltning  dynamical 
systems.  They  are  developing  a  window  interface  that  can  be  used  to  control  computations 
on  a  work  station  or  remote  machine  -  with  the  goal  of  using  powerful  parallel  processors 
as  computational  engines. 

Mahalov  and  Guckenheimer  have  extended  the  work  of  Armbruster,  Guckenheimer 
and  Holmes  on  structurally  stable  homociinic  cycles  to  three-wave  interactions  in  systems 
with  symmetry. 


Berkooz  has  made  rigorous  estimates  using  the  proper  orthogonal  decomposition 
showing  that  a  structured  turbulent  flow,  such  as  the  wall  layer,  has  a  phase  space 
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representation  that  remains  within  a  thin  slab  centered  on  the  most  energetic  modes  for  most 
of  the  time.  However,  exits  from  this  region,  which  is  all  that  our  low-dimensional  models 
include,  .should  not  be  ignored,  since  they  typically  correspond  to  violent  events,  such  as 
the  bursting  phenomenon.  Berkooz  and  Holmes  are  trying  to  develop  a  theory  in  which 
deterministic,  low-dimensional  dynamics  governing  the  low  modes  applies  most  of  the 
time,  passages  from  and  renuns  to  this  being  modeled  probabilistically.  This  might  be 
viewed  as  a  dynamical  closure.  They  plan  to  test  their  theory  on  problems  including  the  32 
and  54  dimensional  projections  of  Aubry  and  Sanghi. 

Sue  Campbell  and  Phil  Holmes  are  continuing  their  studies  of  symmetry  breaking 
(0(2)  — t  D4)  in  systems  with  structurally  stable  heterocliiuc  cycles.  They  have  proved  that 
no  analytic  (second)  integral  of  motion  exists  in  a  certain  limiting  case  and  that  only  two 
pairs  of  the  continuum  of  0(2)  symmetric  heteroclinic  cycles  persist  in  general.  Th^  are 
studying  the  bifurcations  from  these  survivors.  This  work  is  relevant  to  our  nxxiels  of 
interactmg  coherent  structures  in  boundary  layers  with  discrete  spanwise  symmetry,  such 
as  that  caused  by  riblets.  This  is  to  our  knowledge  the  first  analytical  contribution  to  our 
understanding  of  the  drag  reduction  caused  by  riblets. 

Pieter  Swart  and  Phil  Holmes  are  studying  the  dynamic  anti-plane  shear  probkm 
for  nonlinear  viscoelasticity  with  change  (tf  type.  Three  numerical  codes  (finite  difference; 
spectral  and  wavelet  based)  have  been  written  and  are  currently  being  tested.  Theoretical 
studies  of  the  equations,  using  Rybka  &  Kohn's  generalization  of  the  Andrews- Ball-Pego 
transformation,  also  continue.  Although  the  problem  itself  is  only  of  remote  interest  to  die 
subject  ot  this  contract,  the  relative  performance  of  the  three  types  of  codes  is  definitely  of 
interest. 


With  Scott  Jones,  Sid  Leibovich  is  working  on  a  generalized  Lagi^gian  Mean 
decomposition  of  channel  flow.  They  are  first  trying  the  decomposition  on  finite 
amplitude.  2D  Tollmien-Schlichting  waves,  and  then  will  move  to  turbulent  channel  flow. 
The  velocity  fields  are  generated  by  a  DNS  code  developed  by  G.  Yang. 

Leibovich  and  Lumley  proposed  the  application  of  the  averaging  procedure  known 
as  the  Generalized  Lagrangian  Mean  to  the  fully  turbulent  bounda^  layer  problem.  With 
John  Lumley,  Sid  Leibovich  is  developing  a  modelled  set  of  equations  for  ^e  Generalized 
Lagrangian  Mean  of  turbulent  wall  layers.  Leibovich  and  Mahalov  plan  to  extract  cohereut 
structures  from  the  averaged  equations.  The  next  step  would  be  to  study  resonances 
between  these  averaged  structures.  Their  preliminary  estimates  show  that  the  coupling 
coefficients  in  the  equations  describing  interaction  of  average  structures  are  adiabatic 
invariants,  and  are  therefore  important  characteristics  of  the  turbulent  flow.  They  also  hope 
to  be  able  to  show  that  heteroclinic  cycles  are  present  in  the  averaged  (^nations;  these  are 
identified  with  the  bursting  events  in  the  boundary  layer.  Mathematical  work  has  been 
staned  in  this  direction. 

Leibovich  and  Mahalov  have  been  studying  rotating  pipe  flow.  While  this  is  remote 
from  the  subject  of  this  contract,  the  techniques  developed  have  application  to  the  boundary 
layer  (see  below).  Specifically,  they  studi^  interactions  of  distinct  families  of  2-D  helical 
waves,  and  found  a  mechanism,  resonant  triad  interaction,  for  the  interaction  of  these 
waves  resulting  in  fully  three-dimensional  flows. 

Recently  Corke  &  Mangano,  and  Kachanov  have  published  new  dau  on  transition 
in  boundary  layers.  It  appears  that  resonant  triad  interactions  of  Tollmein-Schlichting 
waves  play  a  major  role  here.  Leibovich  and  Mahalov  feel  that  they  will  be  able  to  explain 
many  of  the  results  using  techniques  and  methods  developed  for  rotating  pipe  flow.  One  of 


their  goals  is  to  explain  and  analyze  the  spikes  observed  in  these  expehments.  They  have 
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With  G.  Yang,  Sid  Leibovich  is  anempdng  a  direct  numerical  simulation  of 
turbulence  generated  by  an  oscillating  grid  in  rotating  flows. 
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Research  synopsis 


Current  research  effons  of  Poje  and  Lumley  are  focused  on  accurately  (adequately) 
modelling  the  effects  of  modulated  Reynolds  stresses  on  the  form  of  the  coherent  structures 
predicted  by  an  energy  maximization  principle.  Results  from  use  of  the  anisotropic  eddy 
viscosity  nMdel  show  coherent  structures  which  arc  quite  similar  to  those  predict^  by  the 
Proper  Orthogonal  Decomposition.  Initially,  the  eigenvalue  spectrum  for  these  coherent 
structures  peaked  at  too  low  a  wavenumber,  inclusion  of  feedback  between  the  mean 
velocity  profile  and  the  coherent  structures  moved  the  peak  to  the  observed  'vavenumber. 
We  are  now  extending  this  work  to  include  buoyant  heat  transfer. 

So  far  our  decompositions  have  used  Fourier  modes  in  the  homogeneous  spanwisc 
and  streamwise  directions.  These  distributed  functions,  while  energy-optimal  on  average, 
are  not  well  adaptwi  to  reproduce  spatially  localized  events  such  as  bursts.  In  addition,  one 
of  the  points  on  which  the  work  of  Aubry  et  al  was  criticised,  was  that  the  structures 
produced  were  not  localized  in  space  like  the  structures  observed  in  the  turbulent  boundary 
layer.  Wavelet  bases,  having  compact  suppon,  seem  much  more  appropriate.  J.  Elezgaray 
is  currently  studying,  jointly  with  G.  Berkooz  and  P.  Holmes,  applications  of  wavelets  to 
the  study  of  coherent  structures  arising  from  the  Kuramoto-Sivashinsky  equation.  This 
study  is  preliminary  to  a  more  involv^  study  of  the  more  complex  turbulent  boundary 
layer,  permitting  us  to  gain  experience  with  the  application  of  wavelets  in  chaotic  systems. 
They  have  successfully  modeled  the  wavelet-wavelet  interactions,  and  have  produced  a 
low-dimensional  wavelet  model  that  mimics  in  many  respects  the  behavior  of  the 
Kuramoto-Sivashinsky  equation.  I.  Moroz  and  J.  Lumley  (in  collaboration  with  Berkooz 
and  Elezgaray)  are  working  on  the  extension  of  this  work  to  the  Navier  Stokes  equation. 

To  determine  the  effect  of  control  actuators  on  the  structure  of  the  turbulent 
boundary  layer,  Berkooz,  Holmes  and  Lumley  suggested  a  numerical  scheme  for  a  direct 
(pseudospcctral)  numerical  simulation  scheme.  After  extensive  work  on  the  scheme 
involving  a  linearized  boundary  condition,  it  was  determined  that  it  was  computationally 
unfeasable,  since  the  large  values  of  velocity  near  the  wall  where  the  cell  size  is  small 
resulted  in  unacceptable  time  step  restrictions.  A  new,  exact,  conformal  transformation  (to 
a  flat  boundary)  has  been  devised.  We  have  ported  the  computation  ftom  the  NAS  facility 
at  NASA  Ames  to  the  Cornell  Supcrcomputing  facility;  the  unmodified  code  is  running  on 
both  the  IBM  3090  and  a  cluster  of  RS  6000  machines,  designed  to  exploit  extreme 
parallelism. 

Berkooz  recognized  that  to  get  quantitative  predictions  will  require  generating 
higher  and  higher  dimensional  dynamical  systems,  and  that  the  best  way  to  go  about  this  is 
to  utilize  computer  algebra.  To  this  end,  Berkooz  initiated  an  interdisciplinary  effort 
involving  R.  Zippel  and  S.  Rapkin  from  Computer  Sciences.  The  program  they  have 
written  is  now  capable  of  reproducing  the  work  of  Aubry  et  al  automatically,  carrying  out 
the  Galerkin  projection,  determining  the  coefficients,  and  writing  a  program  for  the  solution 
of  the  ODE'S.  The  work  has  reached  such  a  stage,  that  a  proposal  for  major  support  has 
been  submitted  to  DARPA. 

With  Scott  Jones,  Sid  Leibovich  is  working  on  a  generalized  Lagrangian  Mean 
decomposition  of  channel  flow.  They  arc  first  trying  the  decomposition  on  finite 
amplitude,  2D  Tollmien-Schlichting  waves,  and  then  will  move  to  turbulent  channel  flow. 
The  velocity  fields  are  generated  by  a  DNS  code  developed  by  G.  Yang. 

Leibovich  and  Lumley  proposed  the  application  of  the  averaging  procedure  known 
as  the  Generalized  Lagrangian  Mean  to  the  fully  turbulent  bounda^  layer  problem.  With 
John  Lumley,  Sid  Leibovich  has  developed  a  m^elled  set  of  equations  for  the  GeneralizBd 
[.agrangian  Mean  of  turbulent  wall  layers.  Wing  Ma,  working  with  Leibovich,  has  shown 
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that,  with  a  simple  assumption,  these  equations  will  produce  the  turbulent  mean  velocity 
profile,  as  well  as  the  laminar  one.  This  is  the  first,  lowest  order,  demonstration  that  these 
equations  provide  an  alternate  to  Reynolds  averaging.  Ma  and  Leibovich  are  currently 
using  the  GLM  formulation  to  describe  the  growth  of  Tollmein-Schlichting  waves  in  the 
laminar  boundary  layer.  This  is  a  test  of  the  GLM  formulation,  since  the  results  are  well- 
known  and  well-understood.  When  this  has  been  satisfactorily  resolved,  Ma  will  move  on 
to  the  turbulent  boundary  layer. 

Lumley  considered  the  interaction  between  two  adjacent  bursting  events,  or 
between  an  event  induced  by  an  actuator  (in  an  effort  to  control  the  layer)  and  a  bursting 
event.  He  obtained  expressions  for  the  coefficients  in  the  low-dimensional  model,  as 
nxxlified  by  the  presence  of  the  adjacent  event  The  primary  difference  in  the  coefficients  is 
that  the  coefficient  of  the  linear  term,  representing  pressure  drop  and  loss  to  the  unresolved 
nxxles,  becomes  complex,  causing  the  system  to  spiral  in  toward  a  fixed  point,  instead  of 
falling  directly  toward  it.  Brianno  Coller  and  Phil  Holmes,  starting  from  this  idea,  have 
developed  a  control  algorithm  for  our  simplified  backbone  system,  rotating  the  system 
about  the  fixed  point  first  one  way  and  then  the  other  (or,  easier  to  visualize:  rotating  the 
principal  axes  back  and  forth)  to  keep  the  system  heading  toward  the  fixed  point  as  long  as 
possible.  The  system,  of  course,  finally  escapes  from  control,  and  is  allowed  to  jump. 
This  sort  of  control  results  in  a  very  substantial  delay  in  jumps,  and  an  effective  drag 
reduction  by  something  like  a  factor  of  three.  It  is  somewhat  unrealistic,  however,  in  that 
control  in  a  real  boundary  layer  could  not  be  adjusted  on  an  arbitrarily  shon  time  scale, 
since  one  is  limited  by  the  response  time  of  the  fluid  field.  We  will  look  into  this,  and  also 
develop  a  firm  theOTctical  basis  for  what  is  presently  an  ad-hoc  method. 

Berkooz  has  developtd  a  control  scheme  for  a  low-dimensional  model  of  the 
Kuramoto-Sivashinski  equation,  reminiscent  of  our  low-dimensional  model  for  the 
turbulent  wall  layer.  He  uses  a  control  actuator  which  is  the  analog  of  our  bump  on  the 
wall.  So  far,  the  drag  reduction  achieved  is  of  the  order  of  that  achievable  by  riblets  in  a 
real  boundary  layer.  These  are  preliminary  results,  however,  and  we  believe  that  we  can 
do  better. 
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Lumley,  J.  L.  Four  lectures  at  Universitat  degli  Studi  di  Roma  "La  Sapienza""  Dipardmento 
n.  37  -  Idraulica,  Trasporti  e  Strade:  The  turbulent  boundary  layer.  Drag  Reduction; 
dloherent  Structures;  and  A  low  dimensional  model  of  the  wall  region.  Rome,  Italy,  May 
20-24,  1991.  Invited. 


Lumley,  J.  L.  Low  dimensional  model  of  the  wall  region  in  a  mrbulent  boundary  layer 
new  results,  at  the  lUTAM  Symposium  and  NATO  Advanced  Research  Workshop  on  the 


Interpretation  of  Tune  Series  from  Nonlinear  Mechanical  Systems,  25-30  August,  1991, 
University  of  Warwick,  Coventry,  UK.  Invited. 

Lumley,  J.  L.  Stability,  drag  reduction  and  control  of  the  turbulent  boundary  layer,  using  a 
low-dimensional  model.  Presented  at  the  13th  IMACS  World  Congress  on  Computation 
and  Applied  Mathematics,  July  22-26,  1991,  Trinity  College,  Dublin,  Ireland.  Invited. 

Poje,  A.  &  Lumley,  J.  L.  Coherent  structure  -  fine  grain  turbulence  interaction.  Paper  No. 
BA  1,  Bulletin  of  the  American  Physical  Society  36(10);  2627  .  44th  Annual  Meeting, 
Division  of  Fluid  Dynamics,  American  Physical  Society,  Scottsdale,  AZ,  24-26  November 
1991.  Contributed. 

Ristorcelli,  J.  R.  &  Lumley,  J.  L.  A  rapid  pressure  model  frame  indifferent  in  the  2D  limit 
Paper  No.  HA  6,  Bulletin  of  the  American  Physical  Society  36(10):  2687  .  44th  Annual 
Meeting,  Division  of  Fluid  Dynamics,  American  Physical  Society,  Scottsdale,  AZ,  24-26 
November  1991.  Copntributed. 

Professional  Activiries 

Guckenheimer,  J.  1991.  Editorial  Boards:  Physica  D,  SIAM  Review,  Aequadones 
Mathemadcae,  Intemadonal  Journal  of  Bifurcadon  and  Chaos 

Guckenheimer,  J.  1991.  Director,  Center  for  Applied  Mathematics 

Guckenheimer,  J.  1991.  Directed  of  Research  Programs,  Cornell  Theory  Center 

Guckenheimer,  J.  1991.  Advisory  Board,  The  Geometry  Center  (Minnesota) 

Holmes,  P.  CRM  Special  Year  in  Dynamical  Sytems,  Montreal,  Quebec,  1993-94;  co- 
organizer. 

Holmes,  P.  MSRI  Special  Semester  in  Dynamical  Systems  and  Probabilistic  Method  for 
PDE,  Spring,  1994;  Organizing  committee  member. 

Hoboes,  P.  Editorial  Board  Member.  Archive  for  Rational  Mechanics  and  Analysis; 
Journal  of  Nonlinear  ScinecetNonlinear  Science  Today;  Proc.  Math.  Society  of  Edinburgh. 

Holmes,  P.  Refereeing  Activities;  AMS  Journals,  Applied  Mathematical  Modelling, 
Automadca,  (Communications  in  Mathematical  Physics,  Ergodic  Theory  and  Dynamical 
Systems,  Handbook  of  Acoustics,  Icarus,  Journal  of  Applied  Mechanics,  Journal  of 
Differential  Equations,  Journal  of  Dynamic  Systems  and  Control,  Journal  of  Fluid 
Mechanics,  Journal  of  Mathematical  Biology,  Journal  of  Sound  and  Vibration,  Lecture 
notes  on  Biomathematics,  Mechanics  Research  Communications,  Nonlinearity,  Nonlinear 
Phenomena  (Physica  D),  Physics  of  Fluids,  Physics  Letters,  Rocky  Mountain  Journal  of 
Mathematics,  The  Royal  Society  Journals,  Transpon  Theory  &  Statistical  Physics. 

Holmes,  P.  Reviewing  activities:  American  Scientist,  Applied  Mechanics  Reviews,  Bulletin 
of  the  A.M.S.,  I.M.A.  Journal,  Journal  of  Applied  Mechanics,  Journal  of  Sound  and 
Vibration,  Physics  Today,  Physics  World,  Shock  and  Vibration  Digest,  SIAM  Review. 

Leibovich,  S.  1991.  Associate  Editor,  Journal  of  Fluid  Mechanics 

Leibovich,  S.  1991.  Chairman,  Operations  Committee,  Basic  Engineering  Group 
Operating  Board,  ASME 


Leibovich,  S.  1991.  Chairman,  U.  S.  National  Committee  on  Theoretical  and  Applied 
Mechanics.  (3  year  term). 

Leibovich,  S.  1991.  Co-editor,  ACTA  MECHANICA 

Leibovich,  S.  1991.  Member  of  Executive  Committee  and  Chair  of  External  Affairs 
Subcommittee,  Division  of  Fluid  Dynamics,  Amer  Phys  Soc. 

Leibovich,  S.  1988-1992.  Member,  Congress  Committee  of  the  International  Union  of 
ThMietical  &  Applied  Mechanics. 

Leibovich,  S.  1989-1994.  Member,  Editorial  Board,  Annual  Reviews  of  Fluid  Mechanics 

Leibovich,  S.  1993.  Member,  Scientific  Committee,  2nd  International  Conference  on 
Fluid  Mechanics,  Beijing. 

Leibovich,  S.  1993.  Member,  Scientific  Committee,  lUTAM  Symposium  on  Stability  of 
Nonpaxailel  Flows,  Potsdam. 

Leibovich,  S.  1991.  Member,  Umoshenko  Medal  Selection  Committee,  AS  ME 

Leibovich,  S.  1992.  Member,  United  States  Papers  Selection  Committee,  International 
Congress  of  Theoretical  and  Applied  Mechanics,  Haifa,  ISRAEL. 

Leibovich,  S.  1991.  Member-at-large,  Basic  Engineering  Group  Operating  Board,  ASME 

Leibovich,  S.  1991.  NSF,  Committee  of  Visitors,  Fluid  Dynamics,  Hydraulics,  and 
Particulates  Programs,  Division  of  Chemical  and  Thermal  Systems 

Leibovich,  S.  1991.  Reviewer:  Journal  of  Fluid  Mechanics,  Physics  of  Fluids,  Journal  of 
Physical  Oceanography,  NSF 

Lumley,  J.  L.  Representative  of  Division  of  Buid  Dynamics,  American  Physical  Society  to 
U.  S.  COTgress.  Congressional  Day,  Thursday,  Af^  25. 

Lumley,  J.  L.  Chairman,  Ad-Hoc  Planning  and  Advocacy  Group  for  Turbulence. 

Lumley,  J.  L.  Editor  and  Chairman,  Editorial  Committee,  Subcommittee  on  Research 
Directions  in  Fluid  Mechanics,  U.  S.  National  Committee,  Theoretical  and  Applied 
Mechanics. 

Lumley,  J.  L.  Member,  Scientific  Committee,  lUTAM  Symposium  Eddy  Structure 
Identification  in  Free  Turbulent  Shear  Bows,  Poitiers,  France,  September  14-16,  1992. 

Lumley,  J.  L.  Member,  Advisory  Committee,  8th  Turbulent  Shear  Bows  Symposium. 

Lumley,  J.  L.  Chairman,  Advisory  Committee,  Center  for  Turbulence  Research,  Stanford, 
NASA  Ames. 

Lumley,  J.  L.  Chairman,  Technical  editorial  board,  Izvesnya:  Atmospheric  and  Oceanic 
Physics. 

Lumley,  J.  L.  Co-editor,  Annual  Reviews  of  Fluid  Mechanics. 


Lumlcy,  J.  L.  Editor,  Theoretical  and  Computational  Fluid  Dynamics. 

Lumley,  J.  L.  Member,  editorial  board  Fluid  Mechanics:  Soviet  Research. 

Lumley,  J.  L.  Session  Chairman,  Conference  on  Turbulence  Modeling,  NASA  Lewis, 
August  21,  22. 

Lumley,  J.  L.  Reviewer,  J.  Fluid  Mech.,  Phys.  Fluids,  Physics  Letters  A,  Combustion 
Science  &  Technology,  AIAA  J,  Ira.  J.  of  Multiphase  Flow,  J.  Sci.  Instrum.,  NSF,  ARO 

Lumley,  J.  L.  Expen  wimess.  Breed  Automotive  (Plaintif)  vs.  David  Breed.  Analysis  of 
air  bag  sensor. 

Lumley,  J.  L.  Expert  witness.  Philosopher's  Wool  Co.  (PlaintiO  vs.  Ontario  Hydro. 
Evaluation  of  Hydrogen  Sulfide  monitoring  stations  installed  by  Ontario  Ministry  of  the 
Environment 


